arXiv:1504.05404v2 [cond-mat.mes-hall] 27 Aug 2015 


Tailoring the magneto dynamic properties of nanomagnets using magnetocrystalline 

and shape anisotropies 


Vegard Flovik,^’[^ Ferrari Macia,^ Joan Manel Hernandez,^ Rimantas Brucas,^ Maj Hanson,^ and Erik Wahlstrom^ 

^Department of Physies, Norwegian University of Seienee and Teehnology, N-7491 Trondheim, Norway 
^ Grup de Magnetisme, Dept, de Fisiea Fonamental, Universitat de Bareelona, Spain 
^Department of Engineering Seiences, Uppsala University, SE-751 21 Uppsala, Sweden 
Department of Applied Physies, Chalmers University of Teehnology, SE-412 96 Gothenburg, Sweden 

(Dated: August 31, 2015) 

Magnetodynamical properties of nanomagnets are affected by the demagnetizing fields created 
by the same nanoelements. In addition, magnetocrystalline anisotropy produces an effective held 
that also contributes to the spin dynamics. In this article we show how the dimensions of magnetic 
elements can be used to balance crystalline and shape anisotropies, and that this can be used to 
tailor the magnetodynamic properties. We study ferromagnetic ellipses patterned from a 10 nm 
thick epitaxial Fe him with dimensions ranging from 50 x 150 nm to 150 x 450 nm. The study com¬ 
bines ferromagnetic resonance (FMR) spectroscopy with analytical calculations and micromagnetic 
simulations, and proves that the dynamical properties can be ehectively controlled by changing the 
size of the nanomagnets. We also show how edge defects in the samples inhuence the magnetization 
dynamics. Dynamical edge modes localized along the sample edges are strongly inhuenced by edge 
defects, and this needs to be taken into account in understanding the full FMR spectrum. 


I. INTRODUCTION 

The magneto dynamic properties of nanostructures 
have received extensive attention, from both fundamental 
and applications viewpoint^^^^^. Nanometer sized mag¬ 
netic elements play an important role in advanced mag¬ 
netic storage scheme^^, and their static and most im¬ 
portantly their dynamic magnetic properties are being 
intensely studied. While technological applications are 
important, there is also signihcant interest in understand¬ 
ing the fundamental behavior of magnetic materials when 
they are confined to nanoscale dimensions. In confined 
magnetic elements, there is a complex competition be¬ 
tween exchange, dipolar and anisotropic magnetic ener¬ 
gies. Understanding the interplay between the various 
energy terms is thus of importance when investigating 
the magnetodynamics of such systems. 

The magnetization dynamics in patterned magnetic 
structures has been extensively studied previousl5^^®3. 
The spin dynamics in elliptical permalloy dots were in¬ 
vestigated by Gubbiotti et al^. They studied the various 
excitation modes as a function of dot eccentricity and in¬ 
plane orientation of the applied field, showing how the 
shape of the ellipses affects the spectrum of excitable 
modes and their frequencies. 

However, the above mentioned studies of patterned 
magnetic structures were all performed for systems hav¬ 
ing a negligible magnetocrystalline anisotropy. Mate¬ 
rial systems with a significant crystalline anisotropy pro¬ 
duce an effective field which also contributes to the spin 
dynamics. The combination of shape and crystalline 
anisotropy results in a complex energy landscape, where 
the interplay of these energy terms determines the mag¬ 
netodynamic properties of the system. 

The influence of shape and crystalline anisotropy on 
magnetic hysteresis and domain structures in submicron- 


size Fe particles have previously been investigated by M. 
Hanson et a/.^. However, to the best of our knowledge, 
the dynamic properties of magnetic structures utilizing 
both crystalline and shape anisotropies remains unex¬ 
plored. The goal of this study is thus to investigate a 
system where the energy terms from both crystalline and 
shape anisotropy contribute to determine the dynamics 
of the system. 

We have investigated a system utilizing epitaxial Fe 
as the ferromagnetic (FM) material, patterned to an ar¬ 
ray of elliptical nanomagnets. This results in a system 
combining the cubic crystalline anisotropy of Fe with the 
shape anisotropy due to the elliptical shape of the con¬ 
fined magnetic elements. 

The dynamic properties were investigated by ferromag¬ 
netic resonance (FMR) experiments for ellipses with a 
thickness of 10 nm and lateral dimensions of 50 x 150 
nm, 100 X 300 nm and 150 x 450 nm. The experimen¬ 
tal results are compared with micromagnetic simulations, 
and a macrospin model considering the total free energy 
density of a ferromagnetic structure containing both crys¬ 
talline and shape anisotropies. The macrospin model is 
then used to explore the properties of ellipses with lateral 
dimensions ranging from 50 x 150 nm to 500 x 1500 nm, 
showing how the ellipse size governs the balance between 
crystalline and shape anisotropy. 

During the fabrication of such structures, the magnetic 
properties may be affected by edge defects and shape 
distort ion As the size of the magnetic elements are 

reduced, the edge regions become increasingly important. 
Understanding how edge defects affect the magnetody¬ 
namic properties of the elements is thus of importance in 
nanomagnets, where the edge region covers a significant 
amount of the total sample area. We show how this af¬ 
fects the magnetization dynamics, and that edge defects 
need to be taken into account in understanding the full 
FMR spectrum. 
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II. EXPERIMENTAL SETUP 

The samples are based on a single crystalline Fe film 
epitaxially grown on MgO(OOl) substrates. The ferro¬ 
magnetic ellipses were patterned by e-beam lithography 
and ion beam milling from a 10 nm thick Fe layer, and 
have lateral dimensions of 50 x 150 nm, 100 x 300 nm 
and 150 x 450 nm. The crystalline easy axis [100] and 
[010] of the Fe film are oriented along the long/short axis 
of the ellipses, as indicated in Fig. [TJi. Further details 
concerning sample growth and processing are similar to 
that described earlieJ^. 

The FMR experiments were performed using two com¬ 
plementary setups. The cavity FMR measurements were 
carried out in a commercial X-band electron paramag¬ 
netic resonance (EPR) setup with a fixed microwave fre¬ 
quency of 9.4 GHz (Bruker Bio-spin ELEXSYS 500, with 
a cylindrical TE-011 microwave cavity). The magnitude 
of the external field is then swept to locate the resonance 
field, i^R. The sample is attached to a quartz rod con¬ 
nected to a goniometer, allowing to rotate the sample 360 
degrees in order to accurately resolve the angular depen¬ 
dence. The EMR measurements were performed with a 
low amplitude ac modulation of the static field, which 
allows lock-in detection to be used in order to increase 
the signal to noise ratio. 

Eor the broadband EMR measurements, we used a vec¬ 
tor network analyzer (VNA) EMR setup with a coplanar 
waveguide (CPW) excitation structure. The static exter¬ 
nal field, i^o, was applied in the sample plane, and per¬ 
pendicular to the microwave field from the CPW. This 
was used to obtain the standard microwave S parameters 
as a function of frequency for various fixed values of the 
static field. This allows for a complete field versus fre¬ 
quency map of the resonance absorption, not being lim¬ 
ited to a fixed frequency as for the cavity measurements. 
Data was then collected in a field range of ± 500 mT, 
and a frequency range of 1-25 GHz. Typical absorption 
maps had a step size of A/ = 0.1 GHz and AHq = 5mT. 

III. MICROMAGNETIC SIMULATIONS 

The micromagnetic calculations were performed using 
MuMa:>P^. The simulated ellipses have a dimension of 
150 X 450 nm, with a thickness of 10 nm. In order to have 
mesh independence, the discretization cells should have 
sides of the same order, or less than, the two character¬ 
istic magnetic length scales of the system. The exchange 
length, /exch = magnetostatic exchange 

length Idem = Here A is the exchange stiffness 

constant, Ki the first order anisotropy constant, Kd the 
energy density of the stray field, and an upper limit for 
Kd is given by . 

Material parameters used in the simulations are stan¬ 
dard literature values, with a saturation magnetization, 
Ms = 1.7 X 10^ A/m and a crystalline anisotropy con¬ 
stant of Ki = 4.3 X lO^J/m^, with the easy axis oriented 


along the long and short axis of the ellipse. The exchange 
stiffness was set to a value off A = 21 x 10“^^ Jm“^, and 
the damping coefficient to a = 0.01. 

Performing simulations for a 3d model and a 2d model 
we obtained the same results, and varying the grid size 
it was found that the results converge at a grid size of 
2x2 nm. To save computation time the simulation model 
was thus implemented as a 2d model with a grid size of 
2x2 nm, which is well below the characteristic magnetic 
length scales of the system (/exch = 21 nm. Idem = 3.5 
nm.) 

Simulations of the EMR spectrums were performed by 
using a field relaxation process. The system is first initial¬ 
ized at zero applied field. If a static field, Hq, is applied, 
the simulations are run until the system reaches the new 
ground state configuration. A 10 mT perturbation field, 
Hp is then applied along the z-axis (out of plane), and 
the simulation is run until it reaches the ground state 
configuration for the field Hq + Hp. The perturbation 
field is then switched off, allowing the system to relax. 
The perturbation causes oscillations of the magnetiza¬ 
tion around the equilibrium position with a maximum 
deviation of approx. 1 degrees, avoiding any non-linear 
effects. To obtain the resonance frequencies, we take the 
Eourier power spectrum of the rriz component the first 
10 ns of the magnetization relaxation. The various exci¬ 
tation modes of the system will then appear as distinct 
peaks in the Eourier spectrum.!^ 

Simulations with an ac field of varying frequency as the 
perturbing field were also performed, and we obtained the 
same results as for the field relaxation procedure. The ac 
approach is however more time consuming, as one has to 
scan the full frequency range for each value of the applied 
static field in order to locate all the resonances. To obtain 
the full field vs. frequency map of the excitation modes 
in the system we thus used the field relaxation process. 


IV. FREE ENERGY DENSITY AND 
THEORETICAL FMR SPECTRUM 

Due to the size and shape of the ellipses, we consider 
the individual magnetic elements to be in a single do¬ 
main state. This was also confirmed by MFM imaging of 
similar sample^^, where all particles were found to be in 
a single domain state for a thickness of 10 nm. Increas¬ 
ing the thickness makes it energetically favorable to form 
flux closure domains, and already at a thickness of 30 nm 
some of the particles were found to be in such multi do¬ 
main states. This means that to make sure the magnetic 
elements are in single domain states, one has to keep the 
film thickness well below 30 nm for ellipses of the di¬ 
mensions we have investigated. Having a single domain 
state allows us to use an analytical macrospin model to 
investigate the ferromagnetic resonance properties of the 
system. 

The array of ellipses has an inter-particle spacing of 
two times the corresponding ellipse dimension in each 
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direction, as illustrated in Fig. [T}i. This spacing is suffi¬ 
cient to significantly reduce the dipolar coupling between 
the individual elements, and as a first approximation we 
consider the ellipses as uncoupled magnetic elements. 

We start by defining the geometry of the system, and 
consider the free energy density of the individual mag¬ 
netic elements. From the sample geometry illustrated in 
Fig. (magnetic element in the x-y plane), one gets 
that: 



FIG. 1. a) Array of ellipses with dimension A x an aspect 
ratio of A/B = 3 and inter-particle spacing of two times the corre¬ 
sponding ellipse dimension in each direction. The [100] and [010] 
crystallograpic axis of Fe is oriented along the long/short ellipse 
axis, b) Field geometry of the individual ellipses 


Mx = Ms sin 0 cos 0 

= Mg sin ^ sin 0 ( 1 ) 

Mz = Ms cos^, 

where Ms is the saturation magnetization. Assuming 
the external applied field, Hq, is oriented in the sample 
plane. Oh = 7 r/ 2 , gives 

Hx = Hq cos 
Hy = Hq sin (^H 

After defining the geometry, one can calculate the free 
energy density of the system by adding up the vari¬ 
ous energy terms. Using a macrospin model, we do 
not consider the exchange energy. The total free en¬ 
ergy density of the system is then given by, F^tot = 

A'Zeeman “1“ A'P)g] 2 iagnetization A'^nisotropy • 


Ez = -M -H 

= —MsHq [sin 0 cos (j) cos (pH + sin 0 sin (p sin pn] (3) 
= —MsHq sin^cos (0 — Ph)^ 


^Demag = y + NyMy^ + 




2 , 


Nx sin^ 0 cos^ p E Ny sin^ 0 sin^ p (4) 


-h Nz cos^ 0 


where /ip is the vacuum permeability, Ni are the de¬ 
magnetization factors and Nx + Ny Nz = 1. The units 
for the saturation magnetization and magnetic field are 
[Ms] =A/m and [H] =T respectively. We assume cubic 
crystalline anisotropy for the epitaxial Fe film, with the 
easy axis oriented parallel to the long/short axis of the 
ellipse, as indicated in Fig. The lowest order term in 
the crystalline anisotropy energy is then the fourth order 
term: 


F/Anis — Hl[o(‘^(y^ -|- “h 


= Ki 


sin^ 0 sin^ p cos^ p + sin^ 0 sin^ p cos^ 0 (^ 5 ^ 


+ sin^ 0 cos^ p cos^ 0 


where Ki is the magnetocrystalline anisotropy con¬ 
stant and ai = Mi/Ms. After adding the terms, one 
can write the total free energy density as: 


Etot = -MsHq sin 0 cos{p - pn) 


+ 


JIqMs 


2 

.2 2 , 




sin 0 cos p [ Nx -\ -^ sin 0 sin 


+ sin 0 sin p [ Ny -\ -cos 9 


l^oM^ 
2Ki 2 ^ 


2K 

+ cos^ 0 ( Nz - ^ sin^ 0 cos^ p 

/ioM 2 


(6) 


Equation describes a complex energy landscape, 
with competing energies from the various terms. It is 
important to note that the orientation of the magneti¬ 
zation, given by 0 , might not be parallel to the applied 
field, pH‘ Thus, to investigate the resonance conditions 
of the system one must first find the equilibrium orien¬ 
tation of the magnetization. The equilibrium orientation 
was found by minimizing the free energy density of the 
system given by Eq. for each value of H and pn^ and 
was performed numerically. After obtaining the equilib¬ 
rium orientation of the magnetization, one can calculate 
the resonance frequency cj, given byP^ 


7 ( d^Etot Etot ( d^Etot V\ . 

^ fx^MsSmOX y 80“^ dcfP' \ dOdcf) ) j 

By solving Eq. 0 , one can obtain the resonance fre¬ 
quency as a function of magnitude and direction of the 
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applied field, Calculating the various terms in 

Eq.Q, one obtains: 


d‘^E, 


tot 


a 6»2 


= MgH sin 9 cos{<p — 4 >h) 


4 


cos 20^1 — cos 4^ — 4:fioMgNz/Ki 


2iioM^ 


{N^ + Ny + {N^ - Ny) cos 20)) 


+ (cos 40 + 7) cos 4:0 


(8) 


d^E-, 


tot 


902 


= MsH sin 0 cos(0 — (pn) 


+ 2Ki sin^ 0 cos 40 sin^ 0 


, fioMs\Ny-N,) 

+ --cos 20 


(9) 


d^Etot 

d9d4> 


= MsH cos 0 sin(0 — (j)H) 
+ 8i^i sin 0 cos 0 sin 0 cos 0 
fioM^iNy - N,) 


cos20sin^^ ^]^Q) 


+ 


4i^i 


For thin films, one can simplify these expressions by 
assuming that the magnetization is oriented in the film 
plane, 0 = 7r/2. After introducing the anisotropy field, 
Hk = 2Ki/Ms^ one obtains the resonance frequency 
given by Eq. 0: 


-( 


H cos(0 - 4 >h) + {eIz 

Nx “h Ny -|- (^Nx — Ny^ cos 20 
2 

H cos(0 — (I)h) + Hk cos 40 


)) + f(3 


i40) 


+ iioMg{Ny - Nx) cos 20 


( 11 ) 


Equation (11) gives the resonance frequency for the 


general case, with the assumption that the magnetization 
is oriented in the sample plane. Depending on the shape 
and size of the magnetic elements, one can then adjust 
the demagnetization factors Ni to obtain the resonance 
conditions for various samples. 


In addition to the four-fold symmetry from the cubic 
anisotropy, one notices that in this case there are ad¬ 
ditional terms of two-fold symmetry due to the shape 
anisotropy along the long/short axis of the ellipse. The 
resonance conditions of the system are thus more com¬ 
plicated, and are determined by the interplay of shape 
and crystalline anisotropies. This brings us to the main 
topic of the study, to investigate how tuning the various 
energy terms changes the magnetodynamic properties of 
the system. 


V. RESULTS AND DISCUSSION 
A. Cavity FMR measurements 

The experiments to investigate the angular dependence 
were performed in the X-band cavity FMR setup de¬ 
scribed in section This gives an angular FMR spec¬ 
trum for both the continuous film and an array of ellipses 
of dimension 150 x 450 nm, as shown in Fig. 
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Resonance field, [Tesla], on radial axis Resonance field, [Tesla], on radial axis 
3) vs. rotation in the x-y plane, [degrees] D ) vs. rotation in the x-y plane, [degrees] 


- 0.1 0 0.1 
Field, [Tesla] 


- 0.1 0 0.1 
Field, [Tesla] 


FIG. 2. (Color online) Experimental FMR spectrum for a) con¬ 
tinuous film and b) ellipses of dimension 150 x 450 nm from the 
X-band cavity FMR setup. 


Going from a continuous film to a patterned array of el¬ 
lipses, there is a significant difference. For the continuous 
film, the four-fold symmetry due to the cubic crystalline 
anisotropy in Fe is dominating. For the ellipses the situa¬ 
tion is more complicated, as there are competing energies 
also from the shape anisotropy. 

To investigate this, we compare the experimental and 
theoretical results. By solving Eq.([TT]) after first mini¬ 
mizing the free energy density for each value of H and 
(j)H^ one gets the FMR dispersion relations shown in the 
lower panel of Fig. From Eq.(ll), the relevant pa¬ 
rameters determining the dispersion are the demagneti¬ 
zation factors Ni^ the anisotropy field Hk and the satu¬ 
ration magnetization Mg. In nanometer-dimension mag¬ 
netic structures, estimates of the demagnetization factors 
using an ellipsoidal for mulae are considered to represent 
the anisotropy fields wel P^ * ^^ [ The factors Ni were found 
fronPS, and for an ellipse of dimension 10 x 150 x 450 nm 
they are: Nx ~ 0.005, Ny 0.05 and Nz = 1 — Nx — Ny. 
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The anisotropy field was determined from the ex¬ 
perimental FMR spectrum in Fig. [^, and was found 
to be approx. 50 mT. In the calculations, Mg was ad¬ 
justed to obtain the best fit between the experimental 
and theoretical spectrum, and the best fit was found for 
Ms = 1.5 X 10^ A/m (a reduction of approx. 10% com¬ 
pared to textbook values of Mg for Fe). 

To compare the angular dependence of the theoretical 
spectrum with experimental results from the cavity mea¬ 
surements shown in Fig. one can inverte the solution. 
This rather gives the resonance field Hr, as a function 
of rotation angle for a fixed excitation frequency of 9.4 
GHz, and the inverted solution is shown in the upper 
panel of Fig. To distinguish the effect of crystalline 
anisotropy and shape anisotropy, the same calculations 
were also performed assuming polycrystalline Fe, setting 
i^k = 0. 




FIG. 3. (Color online) Upper figures: Theoretical data for res¬ 
onance field versus rotation angle for a) continuous film and b) 
ellipse of dimension 10 X 150 X 450 nm, with (red) and without 
(blue) crystalline anisotropy. Lower figures: Dispersion for c) con¬ 
tinuous film and d) ellipse of dimension 10 X 150 X 450 nm, with 
(solid lines) and without (dotted lines) crystalline anisotropy. 

Gomparing theory and experiment in Fig. |^and|^ one 
notices that for the continuous film, both show the ex¬ 
pected four-fold cubic symmetry. For the ellipses, the 
theory replicates the ’’heart shape” of the resonance well. 
In the experimental data in Fig. |^, there are also some 
additional weak resonance lines. It is known that regions 
along the sample edges c ould lead to a spectrum of addi¬ 
tional edge modeJ^E^lI^. However, from our experiments 
we observe that the main mode is dominating, and thus 
focus on this in the following. The other resonances are 
characterized and discussed in detail in section IV Cl 


B. Size of the ellipses 

To investigate the interplay of shape anisotropy and 
crystalline anisotropy, we studied ellipses of various lat¬ 
eral dimensions but with the same aspect ratio of 1:3. 
Changing the sample size affects the balance between 
crystalline and shape anisotropy in the free energy den¬ 
sity. As shown using our macrospin model for the main 
FMR mode, this will in turn change the resonance fre¬ 
quency. There are two limiting cases worth noticing: in 
the limit of a very large ellipse, one should expect a be¬ 
havior close to that of a continuous film, where crystalline 
anisotropy is dominating. By gradually reducing the size 
of the ellipse, shape anisotropy becomes increasingly im¬ 
portant. This means that one can use the size of the 
magnetic elements to tune the ratio between crystalline 
and shape anisotropies, and thus change the magnetody¬ 
namic properties of the system. 

Changing the dimensions of the ellipse affects the free 
energy density of the system, given by Eq.(|^. The tran¬ 
sition from a continuous film to a small ellipse can be 
observed by considering the energy landscape of the sys¬ 
tem as a function of the ellipse dimension, as shown in 
Fig.g 

a) Continuous film b) 500x1500 nm Ellipse 
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FIG. 4. (Color online) Free energy density given by Eq. for a) 
continuous film, b) 500 X 1500 nm ellipse, c) 150 X 450 nm ellipse, 
d) 50 X 150 nm ellipse. Film thickness is 10 nm in all cases. 






Figure [^indicates how the free energy density changes 
when one gradually reduces the size of the ellipse from 
the upper limit of a continuous film, to an ellipse of di¬ 
mension 50 X 150 nm. As expected, one notices that in 
all cases the magnetization favors an orientation in the 
sample plane {0 = 90, from sample geometry as defined 
in Fig.f^). For the continuous film and the largest ellipse 
in Fig.]^ and b, one can clearly see the dominating crys¬ 
talline anisotropy, with a four-fold symmetry between the 
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energy minima along the (j) axis. 

In the intermediate case for an ellipse of dimension 
150 X 450 nm, one has two dominating energy minima at 
0 = 0 and 0 = 180 (magnetization along the long axis of 
the ellipse). In addition, there is a quite flat saddle point 
at 0 = 90 (which corresponds to a magnetization along 
the short axis of the ellipse). This is not a stable energy 
minimum, but the flatness of the saddle point means that 
applying a small magnetic field along this axis will create 
a local energy minimum along this direction. 

For the smallest ellipse, the energy landscape is dom¬ 
inated by the two-fold shape anisotropy along the long 
axis of the ellipse. To align the magnetization along the 
short axis of the ellipse (0 = 90) will thus require a quite 
large external field. 

As shown in section |IV| the FMR frequency given by 
Eq.p is determined by the free energy density of the 
system. Adjusting the lateral dimensions of the ellipse 
is thus an important parameter controlling the FMR fre¬ 
quency. From Eq.(pT|, one notices that the resonance 
frequency is determined by contributions of both two¬ 
fold and four-fold symmetry. From this expression, the 
relevant ratio to determine which term will dominate is 
given by H^/iiqMs{N^ — Ny). Changing the ellipse di¬ 
mensions, and thus the demagnetization factors Ni^ af¬ 
fects the resonance frequency significantly, as shown in 
the upper panel of Fig. 

Theoretical data: 



Resonance field, [Tesla], on radial axis 
vs. rotation in the x-y plane, [degrees] 



Experimental data: 

Resonance field, [Tesla], on radial axis vs. rotation in the x-y plane, [degrees] 

C) 150x450 nm 100x300 nm 50x150 nm 



FIG. 5. (Color online) a) Theoretical dispersion for ellipses of 
dimension 50 X 150 nm (Black), 100 X 300 nm (Blue) and 150 X 450 
nm (Red), b) Angular dependence of same data, c) Experimental 
data for ellipse of dimension 150 x 450 nm d) 100 x 300 nm and e) 
50 X 150 nm. 


and 50 X 150 nm. As the dimensions of the ellipse are re¬ 
duced, the two-fold shape anisotropy tends to dominate 
over the crystalline anisotropy, and the ’’heart shape” 
of the spectrum in Fig. IT due to the cubic crystalline 
anisotropy is suppressed. Comparing the theoretical re¬ 
sults with the experimental data in the lower panel of 
Fig-ll they follow the same trend. As the size is re¬ 
duced, the resonance is shifted to slightly higher fields, 
and the ’’heart shape” of the resonance gets suppressed. 

Investigating the opposite limit, one can determine 
when the crystalline anisotropy starts to dominate. Com¬ 
paring the theoretical FMR spectrum for ellipses of di¬ 
mension 150 X 450 nm, 250 x 750 nm and 500 x 1500 nm 
in Fig. one notices that by increasing the size, the ef¬ 
fect of shape anisotropy is suppressed compared to that 
of crystalline anisotropy. 



—500x1500: long axis 
—500x1500: short axis 
— 250x750: long axis 
-250x750: short axis 
- 150x450: long axis 
- 150x450: short axis 
0.2 03 04 0.5 

Applied field, [T] 


b) 

Resonance field, [Tesla], on radial axis 
vs. rotation in the x-y plane, [degrees] 
90 



FIG. 6. (Color online) a) Theoretical dispersion for ellipses of 
dimension 500x1500 nm (Black), 250x750 nm (Blue) and 150x450 
nm (Red), b) Angular dependence of same data. 


For an ellipse of dimension 500 x 1500 nm the dispersion 
starts to look similar along the long/short axis of the 
ellipse, as indicated in Fig. [^. If the only contribution 
was from the crystalline anisotropy, the dispersion should 
be identical along the long/short axis due to the four¬ 
fold symmetry. Comparing the FMR spectrum for the 
largest ellipse in Fig. to that of a continuous film in 
Fig. IT, they look very similar. This indicates that as 
the sample dimensions approach the micrometric scale, 
shape anisotropies play a minor role compared to the 
crystalline anisotropy. 

To summarize the size dependence, we have shown 
that for sample dimensions above approx. 1 /im, crys¬ 
talline anisotropy will dominate. In the opposite size 
limit, shape anisotropy will dominate for sample dimen¬ 
sions below approx. 50 x 150 nm. In this intermediate 
regime, one can thus effectively use the sample size as a 
parameter to tune the balance between crystalline and 
shape anisotropies. 


C. Broadband FMR measurements and 
micromagnetic simulations 


Figure IT and b compare the theoretical FMR spec¬ 
trum for ellipses of dimension 150 x 450 nm, 100 x 300 nm 


The assumption that the magnetization in the individ¬ 
ual ellipses is uniform is a good approximation at the 






















center of the ellipse, but along the edges the magnetiza¬ 
tion will be less uniform due to the demagnetizing fields. 
Regions along the sampl e edg es could lead to a spectrum 
of additional edge modes^^. In addition, there could be 
other spin-wave excitations with non ze ro w ave vectors, 
and correspondingly varying frequencie^^^. To charac¬ 
terize the various resonances, we thus performed a series 
of broadband FMR measurement in combination with 
micromagnetic simulations. 

To obtain a complete field versus frequency map of the 
FMR absorption, we performed experiments using the 
broadband setup described in section [TTj The experimen¬ 
tal FMR absorption peaks were extracted, and are shown 
in the upper panel of Fig.[^ Red dots represent the main 
FMR mode, and the blue squares the additional weaker 
mode. For clarity only a few selected data points are 
included, where the uncertainty in determining the ab¬ 
sorption peak position is of the order of the dot size. The 
experimental results are then compared with the theoret¬ 
ical FMR spectrum from the macrospin model, shown as 
dotted black lines. 
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FIG. 7. (Color online) Upper panel: Experimental FMR spectrum 
for ellipses with field oriented along long/short axis. Experimental 
data is shown as red dots/blue squares, and theoretical spectrum 
as dotted black line, a) Field sweep from negative to positive field 
along long axis, also showing the switching of the magnetization 
at approx. 75mT. b) Similar measurement along the short axis, 
showing the main mode (red dots) and an additional weak reso¬ 
nance at lower frequency (blue squares). Lower panel: Simulated 
FMR spectrum for a single ellipse of dimension 150 X 450 nm with 
the field oriented along the c) long and d) short axis. 


The agreement between theory and experiment is good 
for an applied field oriented along the long axis of the 
ellipse, as indicated in Fig. izl- Sweeping the field from 
negative to positive, one also notices the switching of the 
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magnetization. As the field is swept from negative to 
zero, the FMR frequency decreases as expected. This 
continues also for positive fields until the external field 
is strong enough to overcome the anisotropy favoring the 
magnetization along the long axis of the ellipse. The 
switching is then observed as an abrubt jump in the FMR 
spectrum. 

When applying the field along the short axis there are 
two parallel dispersing lines, as shown in Fig.[^. A high 
frequency resonance and an additional weaker resonance 
at lower frequency, which corresponds well to the addi¬ 
tional resonance also seen in the cavity measurements 
(see Fig. If)- Comparing the measurements along the 
short axis with the theoretical dispersion, one does not 
observe the low field resonance in Fig.[^ (the black dot¬ 
ted line below 100 mT). In this field range the magneti¬ 
zation is not saturated and it is still oriented along the 
long axis of the ellipse, being parallel to the microwave 
(MW) pumping field from the CPW. However, in the 
cavity measurements we observed both resonances, be¬ 
cause the pumping field is, in this case, oriented out of 
the sample plane and thus perpendicular to the magneti¬ 
zation. The first resonance is observed at a field of ^ 50 
mT (see Fig. [^), and a second one at ^ 100 mT, which 
agrees well with the expected resonance fields from the 
theoretical curves shown in Fig. at a frequency of 9.4 
GHz. At higher fields the magnetization in the ellipse 
saturates in the direction of the external field, being per¬ 
pendicular to the MW pumping field from the CPW, and 
thus the theoretical spectrum corresponds well with the 
high-frequency branch of the experimental data. 

Using a macrospin model, one accounts only for the 
main FMR mode. In order to investigate the observed 
low frequency resonance we performed micromagnetic 
simulations. The model was implemented as a single el¬ 
lipse of dimensions 150 x 450 nm with a thickness of 10 
nm, and the simulated FMR spectrums are shown in the 
lower panel of Fig. Comparing the experimental data 
with the micromagnetic simulations, we notice a few dif¬ 
ferences. Applying the field along the long axis of the 
ellipse, the simulated and experimental data both show 
a single dispersing resonance. The simulated FMR fre¬ 
quency is however noticeably higher than the experimen¬ 
tal results. Applying the field along the short axis of the 
ellipse, the differences between the experimental and sim¬ 
ulated FMR spectrum are more significant. The exper¬ 
imental data show two parallel dispersing lines, whereas 
the simulated spectrum shows a whole range of various 
excitation modes. 

A similar splitting of the main mode has been observed 
experimentally in elliptical permalloy dots, and was at¬ 
tributed to a hy bridi zation of the main mode with other 
spin-wave mode^^^. A study of the excitation modes in 
permalloy dots as a function of dot eccentricity has been 
performed by Gubiotti et where they found a large 
range of possible modes depending on the orientation of 
the external field with respect to the axis of the dots. The 
number of modes in our system compared to theirs may 
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be smaller because of the different material parameters 
and sample size. The exchange stiffness in Fe is almost 
twice that of permalloy, and combined with a smaller 
sample size this results in a reduction in the number of 
excitation modes due to the increased exchange energy. 
This was also confirmed in our simulations, where the 
mode splitting disappear when reducing the sample size 
or increasing the exchange stiffness. 

The low frequency branch in Fig. [gi was identified 
by imaging the rriz component from the micromagnetic 
simulations (out of plane component). From the periodic 
oscillations of the magnetization, we determined the low 
frequency resonance to be localized along the edges of 
the ellipse, as indicated in Fig. [^for an applied field of 
150 mT along the short axis. 



FIG. 8. (Color online) ruz component, showing one oscillation pe¬ 
riod of the edge mode at H= 150 mT, corresponding to a frequency 
of approx. 4 GHz. 


After identifying the excitation modes, one needs to 
consider why there is a significant difference between the 
simulated and experimental FMR spectrum. It is known 
that the fabrication process of nanostructures can lead 
to distortions and defects at the sample edge^^^^^. To 
investigate how this would affect the magnetodynamic 
properties, the effects of edge defects need to be taken 
into account in the simulation model. 


1. Edge modes and edge defects 

In the initial simulations, the edges of the ellipse 
were treated as ideal. However, the samples most likely 
have some kind of non-ideal edges which could influ¬ 
ence the FMR spectrum. The effects of non-ideal edges 
on the dynamics have been investigated theoretically by 
McMichael et al^. It was shown that several cases such 
as edge geometry, reduced edge magnetization and sur¬ 
face anisotropy on the edge surface all had similar effects. 
The main effect was to reduce the edge saturation field, 
which is the field needed to align the magnetization at 
the edge nearly parallel to the applied field. A reduced 
edge magnetization will also lead to a smaller effective de¬ 
magnetization field along the edges. This would cause a 
significant increase of the edge mode resonance frequency 
compared to that of an ideal edge, and the shift could be 
in the order of several ghM Such effects would be less 


important when the field is oriented along the easy axis of 
the ellipse, explaining the better agreement between the 
simulated and experimental spectrum in that geometry. 

To account for edge defects in the simulations, we made 
a model where the material properties were changed 
along the edges of the ellipse. In a real sample the vari¬ 
ation of the material properties when approaching the 
sample edge should be gradual, but as a first approxi¬ 
mation the model was defined with two distinct regions. 
The width of the edge region was set to 10 nm, and is 
within the same width range as that investigated theoret¬ 
ically by McMichael et al?^. A schematic of the model 
including edge defects is shown in Fig. [^. 

As mentioned in section |TI| the samples were defined 
by ion beam milling. This can affect the magnetic prop¬ 
erties of the sample^, and a more disordered edge region 
could lead to an increased damping of the FMR modes. 
Two kinds of defects have thus been considered in the 
simulations; increased damping a, and reduced M^. 



C) Edge: Reduced Ms d) Edge: Increased damping 
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FIG. 9. (Color online) a) Simulation using normal edges, showing 
the edge mode and splitting of the main mode, b) Schematic of 
simulation model with a defined edge region, c) Simulation with 
reduced Mg in edge region (reduction of 40 %). d) Simulation with 
increased damping a in edge region (from 0.01 to 0.1). 


In the initial simulation model with ideal edges, excited 
spin waves would be reflected at the edges of the sample. 
This explains the multiple excitation modes observed in 
the simulations, due to a hy bridi zation of the main mode 
with other spin-wave mode^^^ (see Fig. §1). As a dis¬ 
ordered edge region could lead to increased damping of 
the FMR modes, we introduced an edge region where the 
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damping was increased from a = 0.01 to a = 0.1. This 
would absorb the propagation of spin waves, reducing 
the spin-wave reflection at the sample edges. As seen in 
Fig. IJi, the increased damping lead to a broadening of 
the FMR modes, and suppress some of the splitting of 
the main mode. The low frequency edge mode however, 
remains relatively unaffected. 

The edge magnetization Mg was found to be the most 
important parameter, and we made simulation models 
where the outer region of the ellipse had a significantly 
reduced Mg from Mg = 1.7 x 10^ A/m to Mg = 1 x 10^ 
A/m. Reducing Mg in the edge region changes the FMR 
spectrum considerably, as seen when comparing Fig. 
and c. The splitting of the main mode is suppressed, 
and the resonance frequency of the edge mode is shifted 
significantly. The resulting spectrum now resembles the 
experimental data, showing mainly two parallel dispers¬ 
ing modes. 

Another important effect to consider in arrays of nano¬ 
magnets is the dipolar interaction among the individual 
particles. In order to take this into account, we per¬ 
formed simulations for arrays of interacting ellipses. 


2. Dipolar interactions 

The simulations so far have been performed for single 
ellipses. However, due to the periodic array of ellipses (as 
shown in Fig. Ef), there will be some degree of dipolar 
interaction between the individual ellipses. The dipolar 
interaction in arrays of magnetic particles can have both 
static and dynamic contributions. The effects of static 
dipolar interaction on the magnetization reversal of the 
same samples have been investigated previously, and an 
interaction field in the order of tens of mT was founcP^. 
The dynamic interaction can couple the magnetization 
dynamics of adjacent dots through the stray field gener¬ 
ated by the precessing magnetization, forming collective 
spin excitations in the systenJ^. 

Interactions were included in the simulations by using 
periodic boundary conditions (b.c.), with the same pe¬ 
riodicity as that indicated in Fig. In the limit of 

strong dipolar interaction, one could also expect collec¬ 
tive modes in the system. A simple model of a single 
ellipse with periodic b.c. would not be sufficient to re¬ 
solve such modes, as the neighboring ellipses could rotate 
either in p hase (acoustic mode) or out of phase (optic 
modep^^. To take this into account, we compared the 
simulation results for a single ellipse with periodic b.c. 
versus arrays of 3 x 3, 5x5 and 10 x 10 ellipses. Compar¬ 
ing the simulated FMR spectrums for the various array 
sizes, we found no indication of such collective modes in 
our system. In the following simulations the dipolar in¬ 
teraction was thus taken into account by using a simple 
model for a single ellipse with periodic b.c. 

Comparing the simulated spectrums for a single ellipse 
versus an array of ellipses, we found that the dipolar in¬ 
teraction changes the effective field felt by the individual 


ellipses. At zero applied field, the magnetization is ori¬ 
ented along the long axis of the ellipses. The overall 
dipolar field caused by the array geometry will then op¬ 
pose the magnetization direction. As seen in Fig. 
the dipolar field reduces the resonance frequency at zero 
applied field for the array compared to a single ellipse. 
Increasing the field along the short axis of the ellipse, the 
magnetization will reorient itself along the short axis at 
an applied field of approx. 75 mT (seen as a ’’dip” in 
the FMR spectrum in Fig. [^ ). At fields above this 
switching field, the dipolar interaction acts to increase 
the effective magnetic field felt by the ellipses, and thus 
increases the FMR frequency. These shifts can be seen in 
Fig. for an applied field between 150 mT - 350 mT, 
and are in the order of 1 GHz. These shifts in the FMR 
frequencies along the hard/easy axis are similar to those 
observed by Carlotti et who studied the effects of 

dipolar interactions in arrays of rectangular permalloy 
dots. 

To capture all significant effects we thus made a simu¬ 
lation model with periodic b.c., where edge defects were 
modelled as a reduced Mg at the sample edges. After 
including both edge defects and dipolar interactions, one 
can compare the simulated and experimental spectrums 
in Fig. [Top and c. 
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FIG. 10. (Color online ) a) FFT spectrum of a single ellipse (black 
dots) vs. an array of ellipses (blue line) using periodic b.c, for an 
applied field oriented along the short axis of the ellipse, b) Left 
half: Experimental data and macrospin model as dotted black line. 
Right half: Micromagnetic simulation including edge defects and 
dipolar interactions in an array of ellipses. Data shown for field 
oriented along the long axis, c) same data for field oriented along 
the short axis 
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We notice that the inclusion of edge defects and dipolar 
interactions give a better agreement between the simu¬ 
lated and experimental FMR spectrum. As expected, 
the edge modes are strongly influenced by edge defects 
in the samples. To accurately capture the behavior of all 
the FMR modes, it is thus important to take edge defects 
into account in the simulation model. Due to the large 
spacing between the individual ellipses in the array, the 
dipolar interaction is quite weak. In the simulations we 
observe a small shift in the FMR frequencies, but not 
any indications of collective modes between neighboring 
ellipses. 

The fact that the amplitude of the main mode domi¬ 
nates in the experiments, together with the weak dipolar 
coupling, explains the good agreement between the an¬ 
alytical macrospin model and experimental data. This 
indicates that in the limit of weak dipolar interaction, 
our macrospin model can be used to estimate the FMR 
frequency of the main mode in magnetic elements within 
the investigated size range (e.g. in a single domain state). 
Using an analytical macrospin model compared to per¬ 
forming numerical simulations simplifies the analysis con¬ 
siderably. The various energy terms contributing to the 
FMR dynamics can then be separated, and their relative 
importance investigated. 

VI. CONCLUSIONS 

In this study we have investigated how the com¬ 
bined interplay between shape anisotropy and crystalline 
anisotropy affects the magnetodynamic properties of con¬ 
fined magnetic elements. We have shown how the dimen¬ 
sions of the magnetic elements can be used to balance 
crystalline and shape anisotropies, and that this can be 
used to tailor the magnetodynamic properties 

We have shown that a simple macrospin model for 
the FMR frequency gives good agreement with the ex¬ 
perimental results for the main FMR mode. Com¬ 
paring experimental data and model calculations, we 
show how changing the sample size affects the magne¬ 
todynamic properties. For the smallest ellipses, shape 
anisotropy is dominating, whereas for the largest ellipses 


crystalline anisotropy is the dominating energy term. 
From Eq.(pT]), the relative contributions to the resonance 
frequency from crystalline and shape anisotropy is given 
by: Hk/jiioMs{Nx — Ny), determined by the anisotropy 
field the saturation magnetization Ms and the de¬ 
magnetization factors Ni. This means that for the case 
of a 10 nm thick epitaxial Fe film, one has an intermedi¬ 
ate regime between approximately 50 nm to 1 jam where 
one can use the sample size as an additional tuning pa¬ 
rameter for the dynamic properties. For other materials 
with a different and M^, this regime can be shifted 
to smaller/larger sample sizes. 

The effects of non ideal sample edges and dipolar in¬ 
teraction in the array of ellipses were investigated using 
micromagnetic simulations. We found that edge defects 
in the form of a reduced edge magnetization had to be in¬ 
cluded in the micromagnetic model, and that this needs 
to be taken into account in understanding the full FMR 
spectrum. The static dipolar interaction in the array was 
found to shift the FMR frequency in the order of 1 GHz 
compared to that of a single ellipse. From the simulated 
FMR spectrums we found no indications of collective spin 
excitations due to the dynamic dipolar interaction be¬ 
tween neighboring ellipses. 

The tunability of the relative contributions from crys¬ 
talline and shape anisotropies means that by changing 
the material parameters and sample size one can tai¬ 
lor the magneto dynamic properties of the magnetic el¬ 
ements, which could be of importance for magnonics ap¬ 
plications. 
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